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Possible mechanisms of electronic phase separation in oxide interfaces 
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LaA10s/SrTi03 and LaTi03/SrTi03 (LXO / STO) interfaces are known to host a strongly inho¬ 
mogeneous (nearly) two-dimensional electron gas (2DEG). In this work we present three unconven¬ 
tional electronic mechanisms of electronic phase separation (EPS) in a 2DEG as a possible source 
of inhomogeneity in oxide interfaces. Common to all three mechanisms is the dependence of some 
(interaction) potential on the 2DEGs density. We first consider a mechanism resulting from a siz¬ 
able density-dependent Rashba spin-orbit coupling. Next, we point that an EPS may also occur 
in the case of a density-dependent superconducting pairing interaction. Finally, we show that the 
confinement of the 2DEG to interface by a density-dependent, self-consistent electrostatic potential 
can by itself cause an EPS. 


I. INTRODUCTION 

The observation of a two-dimensional (2D) metallic 
state at the interface of two insulating oxides^, and 
the subsequent demonstration of its gate-tunable metal- 
to-superconductor transition 3 ®, have attracted much at¬ 
tention in the last decade. Numerous experi ments , like 
transport^®], magnetometr}ffin3 tunnelin g 16 * 17 *, and 
piezo-force spectroscopy^, indicate that the 2DEG is 
inhomogeneous. It seem s likely that inhomogeneities 
at nanometric scaled^^D coexist with structural inho¬ 
mogeneities at micrometric scales 20 . Although extrinsic 
sources like local or extended defects may contribute to 
this inhomogeneity, the dense nanoscopic character of the 
electronic charge clearly indicates that an intrinsic mech¬ 
anism is at work leading to a nanoscopic EPS and to 
an inhomogeneous density distribution of the electrons. 
In the present paper, we give an overview of three such 
mechanisms which could account for these generic fea¬ 
tures of the oxide interfaces. Based on the current state of 
knowledge it is hard to gauge the individual relevance of 
each mechanism (although they can work cooperatively). 
Therefore, we deem a systematic investigation worth be¬ 
ing pursued. The paper is structured as follows: In Sec. 
2 we present a mechanism based on the density depen¬ 
dence of the Rashba spin-orbit coupling (RSOC), which 
was found to be sizable in these systems. The second 
mechanism is based on the occurrence of superconduc¬ 
tivity upon varying the electron density and is described 
in Sec. 3. In Sec. 4 we present the third mechanism 
based on the density dependence of the electrostatic po¬ 
tential confining the electrons at the interface. Our final 
considerations are contained in Sec. 5. 


II. ELECTRONIC PHASE SEPARATION FROM 
RASHBA SPIN-ORBIT COUPLING 


We consider a 2d lattice model with strong RSOC de¬ 
scribed by the hamiltonian 
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Here the first term describes the kinetic energy of elec¬ 
trons on a square lattice with nearest-neighbor ((ij)) hop¬ 
ping. The second term is the RSOC 

H RSO __ [7 i,i+yji,i+y - li,i+xjli +x ] (2) 


where = -i'E aa , [4a T aa' c i+v^ ~ 4+^^'^ 

denotes the ^-component of the spin-current flowing on 
the bond between Ri and Ri+rj. 

Following Ref. 2 ^we assume that the coupling constants 
depend on a perpendicular electric field E which is pro¬ 
portional to the local charge density. Since in real space 
the coupling constants 7are defined on the bonds, we 
discretize E at the midpoints of the bonds and define the 
dependence on the charge as E i+r} / 2 = eo + ei(n* + n^ +?7 ). 
For the dependence of the RSOC on the electric field we 
adopt the form given in Ref. 21 ^so that altogether the fol¬ 
lowing coupling is considered 
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For strong RSOC this coupling will induce the forma¬ 
tion of electronic inhomogeneities and thus concomitant 
variations in the local chemical potential \ obtained self- 
consistently by minimizing the energy with respect to the 
density: 
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The left panel of Fig. [l] reports an example of EPS. A 
small amount of disorder is introduced to make the con¬ 
vergence to the inhomogeneous state easier, but in the 
absence of a density dependence of the RSOC it would 
have a minor effect leaving the system essentially homo¬ 
geneous. In this latter case the ground state is character¬ 
ized by a homogeneous flow of j x ^ spin currents along 
the y-(x)-direction whereas the z-component j z vanishes. 
On the other hand, the presence of a density dependent 
RSOC gives rise to a strongly inhomogeneous state which 
induces a finite flow of j z spin currents (black arrows). 
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FIG. 1. Left panel: Charge distribution (red circles) and 
spin current (z-component) in the ground state. Disorder 
potential Vo = 0.It, n e i = 26 particles on a 16 x 16 lattice 
with periodic boundaries. The RSOC is from Eq. § with 
ao = 0.5, a\ — 1.5, fio = /3\ — 0. Right panel: Charge 
distribution (red circles) and current (arrows) in the ground 
state for N = 22 electrons at B — 1692T and open boundary 
conditions with ao — 0.5, a\ — 1.5, /3o = /5i = 0. 


These are caused by the in-plane electric fields E, y=x ,y 
due to the inhomogeneous charge distribution which ac¬ 
cording to the response equation jp = a s e OL ^ 1 E 1 (s a p 7 
denotes the totally antisymmetric tensor and cr s the spin- 
Hall conductivity) generate j z spin currents along the 
x— and y— directions. In general these currents are not 
conserved because the inhomogeneous RSOC introduces 
torque terms which act as source or drains of the spin 
current 23 . Analogous non-conserved spin-currents bound 
to the inhomogeneous charge density are observed for j x 
and j y (not shown). The right panel of Fig. [l] displays 
an inhomogeneous state for the same model of Eq. (lj 
in the presence of a strong perpendicular magnetic field. 
The system is separated in a phase of density ~ 0.31 
and the vacuum phase {n\ = 0). Quite interestingly we 
found that in the quantum Hall regime the EPS instabil¬ 
ity is realized even for very small RSOCs. The reason for 
this surprising result is that the EPS more easily occur in 
states with high density of states (DOS). This naturally 
favors the instability for the Landau-Hofstadter states, 
which have a delta-like DOS. The fact that the (density- 
dependent) RSOC produces an inhomogeneous state in 
the QH regime is an intriguing result that naturally calls 
for the investigation of edge states around and inside the 
inhomogeneous stated. 


III. ELECTRONIC PHASE SEPARATION 
FROM DENSITY DEPENDENT 
SUPERCONDUCTIVITY 

A. Phenomenological Ginzburg-Landau model 

A very interesting feature of LXO/STO interfaces is 
that superconductivity may be tuned by a gating poten¬ 
tial which varies the density of the 2DEG. Therefore these 
systems are characterized by a density dependent super¬ 
conducting critical temperature T c (n). We consider a 



FIG. 2. (a) Chemical potential vs. density for various values 
of the parameter r) in T c (n) — rj(n—n c ) and n c = 0.026 for the 
GL model; (b) Chemical potential and SCOP as a function of 
the electron density n, with n c = 0.04, a — 0.1 and go — 2.7 1 
for the model of Eq. <§• 


standard Ginzburg-Landau (GL) model in which the su¬ 
perconducting order parameter (SCOP) is coupled to the 
density via the density dependence of the critical temper¬ 
ature. Assuming a spatially constant A one has 

F [A] = a \T — T c (n)\ A 2 + 6A 4 . (5) 

Considering the simple form for T c (n) reported in the 
inset of Fig. [2] (a), which starts linearly at n = n c , 
T c (n) = r](n—nc ), and then saturates to a constant value, 
one finds for the SCOP 

A2(n,T) = ^ (Tc(n) _ • 

In the absence of SC the chemical potential is a smoothly 
increasing function (positive electron compressibility) 
that can be obtained, e.g., from a generic van der Waals 
model for the 2DEG. Taking the standard expression of 
the chemical potential fi = dF/dn one obtains that this 
smooth chemical potential is modified once n > n c . Fig. 
|2ja) reports /i vs. n for various values of the parameter 
rj. Clearly, if rj > r\ c fi{n) acquires a negative slope, i.e. a 
negative compressibility, marking a PS instability region, 
which must be found with a standard Maxwell construc¬ 
tion. Clearly the PS dome shrinks upon decreasing r] and 
vanishes at n — n c — 0.026 


B. Microscopic model with density dependent 
pairing 

The phenomenological Ginzburg-Landau approach can 
be complemented with a model, where electrons, de¬ 
scribed by a tight-binding model (for simplicity we con¬ 
sider only a nearest-neighbor hopping t on a square lat¬ 
tice), feel a pairing potential g(n). A simple phenomeno¬ 
logical form which reproduces a rapid growth of a SCOP 
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above a critical density can be introduced, e.g., as 


g(n) = go tanh 


n — n c 
an 


x 6{n — n c ). 


(6) 


go is a parameter controlling the maximum value attained 
by g after its saturation and a is a constant factor which 
determines how steep is the slope of g(n) [and conse¬ 
quently of A(n)]. At any fixed average density one can 
solve the coupled equations determining the SCOP (the 
standard BCS equation) and the chemical potential. Fig. 
[ 2 ] (b) reports the behavior of the SCOP and of the chem¬ 
ical potential g as a function of the electron density near 
the critical density n c = 0.04 for go = 2.7 1. Interestingly, 
the chemical potential displays a non monotonic behav¬ 
ior, marking a region with negative compressibility start¬ 
ing in the proximity of the rapid growth of the SCOP. Not 
only this finding agrees with the scenario emerging from 
the phenomenological GL approach, but one finds that 
the EPS already occurs in a region of relatively weak 
superconducting coupling. Indeed, since the electronic 
tight-binding band has a DOS N(Ep) ~ l/8t, the BCS 
coupling corresponding to go = 2.7 1 is in a regime of 
weak/moderate coupling A ~ 0.34 at which the chemical 
potential barely starts to depend on A. This finding is 
interesting because it indicates that the EPS mediated 
by a density-dependent pairing does not require strongly 
coupled electron pairs. 


IV. ELECTRONIC PHASE SEPARATION 
FROM INTERFACE CONFINEMENT 

Another well known mechanism providing a negative 
contribution to the electronic compressibility is given by 
the mechanism confining the EG at interfaces. In this 
latter case, the self-consistent solution of the Schrodinger 
and Poisson equations relating the electronic wavefunc- 
tions and the electric potential arising from external po¬ 
tential and electronic density itself usually shows that 
the EG becomes more compressible once its finite trans¬ 
verse confinement is taken into account. The large di¬ 
electric constant of the material hosting the EG favors 
its softening, leading to an increased compressibility. In 
this regard, the STO, having a huge dielectric constant 
e > 300 is an optimal candidate to investigate whether 
this mechanism may lead by itself to an overall negative 
compressibility and consequently to an EPS. Indeed, by 
numerically solving the coupled Schrodinger and Poisson 
equations to determine the electronic density and sub¬ 
band structure, we find that this is the case. Fig. [ 3 ] (a) 
schematically displays the profile of the potential well 
along 2 perpendicular to the LXO/STO interface. Upon 
increasing the density, an increasing number of levels is 
filled, but the self-consistent electronic structure adjusts 
to lower its energy and, above some density, it displays a 
downward drift [Fig. ^b)} and a broadening when sub¬ 
bands involving more z-oriented orbitals (like the d xz and 
d yz of Ti) [Figj3|c)]. These self-consistent modifications 



FIG. 3. (a) Schematic behavior of the selfconsistent electro¬ 
static confining potential and of the sub-band electronic states 
perpendicularly to the LXO/STO interface. The LXO layer 
contains the positive countercharges (coming from the polar¬ 
ity catastrophe and/or from oxygen vacancies), while the oc¬ 
cupied sub-band levels are in marked in green, (b) When the 
density of electrons and countercharges is (locally) higher the 
potential well becomes deeper (b) and when d XZjyz orbitals of 
Ti in STO start to be filled (c) the well also becomes broader. 
This induces a decrease of the chemical potential g (dashed 
red line). 



FIG. 4. Density dependence of the Fermi energy (solid blue 
curve), of the chemical potential (solid red curve), of the 
d xy levels (dashed grey curves), and of the d x z or d yz lev¬ 
els (dashed black curves). 


of the potential and of the electronic levels entail a pulling 
down of the chemical potential which eventually gives rise 
to a negative compressibility. Fig. [4] displays the elec¬ 
tron density dependence of the sub-band Ti t^g levels, 
of the Fermi level, and of the chemical potential. Notice 
that the compressibility stays negative [negative slope of 
g(n)\ up to unphysically large densities. Of course several 
physical mechanisms will eventually stop this unphysical 
charge segregation (see below). 


V. DISCUSSION AND CONCLUSIONS 

All the models outlined above are based on the com¬ 
mon idea that some interactions lead to an energy de¬ 
crease when the electron density is increased. In the pres¬ 
ence of a RSOC 7 the bottom of the electronic band is 
lowered by a quantity eo ~ y 2 . Thus, if 7 increases with n 
as in Eq. Q, there is an energetic benefit for the system 
to increase (locally) the density by segregating electrons. 
The same occurs in the model with density dependent 
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pairing §): If the electronic density is increased above 
the threshold n c , electrons can be paired gaining bind¬ 
ing (and possibly condensation) energy. In the electro¬ 
static confinement model an energetic gain occurs when 
the electrons and the positive countercharges in the LXO 
layer (those due to the polarity catastrophe and/or oxy¬ 
gen vacancies) are denser. In this case (cf. Fig. [3] (b)) the 
electrons at the interface feel a deeper confining potential 
and lower their energy by phase separating. 

All the above mechanisms of EPS are obviously con¬ 
trasted by the Coulombic force that opposes the segre¬ 
gation of charged electrons. However, two possibilities 
remain open to allow for an inhomogeneous electronic 


distribution. On the one hand the electrons may sepa¬ 
rate on finite scales until the energetic gain is compen¬ 
sated by the cost of the Coulombic repulsion. Simple 
estimated^ show that the very large value of the dielec¬ 
tric constant of STO weakens the Coulomb repulsion and 
allows this frustrated EPS mechanism to produce rather 
large (^50 nm) inhomogeneities. On the other hand, 
it is also possible that the positive countercharges (like 
the oxygen vacancies) diffuse and follow the segregating 
electrons keeping charge neutrality. Of course also in this 
case EPS stops when the segregating electrons become 
too dense for the countercharges to follow, but finite in¬ 
homogeneities of substantial size can still be formed. 
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